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C~^ • Abstract. We apply the Constitution compilation of 397 supernova la, the baryon 

^^ , acoustic oscillation measurements including the A parameter, the distance ratio 

and the radial data, the five-year Wilkinson microwave anisotropy probe and the 

Hubble parameter data to study the geometry of the universe and the property 

k><( , of dark energy by using the popular Chevallier-Polarski-Linder and Jassal-Bagla- 

^ ' Padmanabhan parameterizations. We compare the simple x^ method of joined contour 



estimation and the Monte Carlo Markov chain method, and find that it is necessary 
to make the marginalized analysis on the error estimation. The probabilities of 
ilk and Wa in the Chevallier-Polarski-Linder model are skew distributions, and the 
marginalized la errors are fl^ = 0.279l°;J5J^, Qk = 0.005to;™^, wq = -1.05+j;;g3, 
and Wa = 0.5l5'5- ^*^^ ^^^ Jassal-Bagla-Padmanabhan model, the marginalized la 
errors are flm = 0.28ltom^, ^k = O-OOOto^aol^ wq = -OMtoil, and Wa = -O-^til- 
The equation of state parameter w{z) of dark energy is negative in the redshift range 
< 2; < 2 at more than 3a level. The flat ACDM model is consistent with the current 
observational data at the la level. 



PACS numbers: 98.80.-k,98.80.Es 
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1. Introduction 

The accelerating expansion of the universe was first discovered by the type la supernova 
(SN la) observations |1] |2]. The phenomena of acceleration could be explained 
straightforwardly by introducing an exotic source of matter with negative pressure, 
the so-called dark energy, which dominates the total matter content of the universe at 
the present epoch and causes the expansion to accelerate. During the past decade, in 
addition to the simple cosmological constant model, a lot of dynamical dark energy 
models, such as the quintessence [3], phantom [3], k-essence [5], tachyon |6], quintom 
[7J, h-essence p], Chaplygin gas [9], holographic dark energy [10], f{R) [TTj, Dvali- 
Gabadadze-Porrati |l12j models, etc, have been proposed. Although a lot of efforts 
have been made to understand the driving force of the accelerating expansion and 
the property of dark energy, whether dark energy is dynamical or not is still an open 
question. Therefore, it is necessary to study the nature of dark energy such as the 
evolutions of its energy density and equation of state. 

Apart from phenomenological models, another effective approach to study dark 
energy is through the observational data. Recently, based on the popular Chevallier- 
Polarski-Linder (GPL) parametrization of dark energy [13], it was found that the fiat 
AGDM model is inconsistent with the current data at more than la level [HI [15]. In 
[14j, it was suggested that the cosmic acceleration is slowing down from z ~ 0.3. In [T5] . 
it was claimed that dark energy suddenly emerged at redshift z ~ 0.3. Furthermore, 
possible oscillating behavior of dark energy was found in [16] . However, no evidence for 
dark energy dynamics was found in [T7| [T8| [19] . It was argued that the systematics in 
different data sets heavily affected the fitting results from observational data [TH [19] . To 
further study the dynamics of dark energy, it is necessary to apply more complimentary 
observational data. In this paper, we combine the Gonstitution sample of 397 SN la data 
[20j, the model independent A parameter from the baryon acoustic oscillation (BAG) 
measurements [21], the two BAG distance ratios at z = 0.2 and z = 0.35 [22], the radial 
BAG measurements aX z = 0.24 and z = 0.43 [23], the five-year Wilkinson microwave 
anisotropy probe data (WMAP5) [23], and the Hubble parameter H{z) data [25l[26] to 
probe the geometry of the universe and the nature of dark energy by using the GPL and 
Jassal-Bagla-Padmanabhan (JBP) [27] parameterizations. We first use the simple x^ 
method of joined contour estimation to obtain the constraints on the model parameters. 
However, the simple x^ method by fixing other parameters at their best fit values has 
some drawbacks because we neglect the correlation effects between the parameters and 
the degeneracy between parameters was not considered. When the parameters are 
strongly correlated, the errors of some parameters will be under-estimated if we fix the 
other parameters at their best fit values. So we also apply the Monte Garlo Markov 
chain (MGMG) method to obtain the marginalized errors of the model parameters. The 
advantage of the MGMG method is that it considers the correlations between the model 
parameters and the result is more reliable. 

The paper is organized as follows. In section 2, we present the SN la data [20], the 
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BAO data [211 El E3] , the WMAP5 data [21] and the H{z) data, and all the formulas 
related with these data. In section 3, We use the ACDM model as an example to show 
how to apply the data to constrain cosmological models. In section 4, we use the CPL 
model to study the geometry of the universe and the property of dark energy. The JBP 
model is used to probe the geometry of the universe and the evolution of dark energy 
in section 5. We conclude the paper in section 6. 



2. Fitting procedure 

To use the Constitution compilation of 397 SN la data [20], we minimize 



2 ^ Y^lMfi) - ^^ohs{zi)f 
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where the extinction-corrected distance modulus /i(z) = 5\ogiQ[dL{z) /Mpc] + 25, cTi is 
the total uncertainty in the SN la observation, the luminosity distance dii^z) is 
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the dimensionless Hubble parameter E{z) = H{z)/Hq, and 
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Due to the arbitrary normalization of the luminosity distance, the nuisance parameter h 
in the SN la data is not the observed Hubble constant. So we marginalize the nuisance 
parameter h with a flat prior, after the marginalization, we get [28], 
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where aj = fiobs{zi) — 25 — 5 logiQlHodi^Zi)], and p denotes the fitting parameters in the 
model. When using the SN la data, the radiation term can be neglected because its 
contribution is negligible. 

In addition to the Constitution SN la data, we use the BAO distance measurements 
from the oscillations in the distribution of galaxies. From the BAO observation of the 
galaxy power spectra, Percival et al measured the distance ratio 

rs{zd) 
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at two redshifts z = 0.2 and z 

0.1097 ± 0.0036, respectively [22]. Here the effective distance is 
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the drag redshift z^ is fitted as 






the comoving sound horizon is 

E(z 



rs{z)=l ^#?, (9) 



the sound speed c,(z) = 1/^3[1 + Rb/{1 + z)], and Rb = 3Qbh^/{A x 2.469 x 10"^). To 
use these BAO data, we calculate 

xIao2(P5 ^6^^ h) = AxiCovi'^{xi, Xj)Axj, (10) 

where Xj = {dz=o.2,dz=o.35), Axj = Xi — xf"^ and Covi(xj,Xj) is the covariance matrix 
for the two parameters (io.2 and (io.35 [22]. Besides the model parameters p, we need to 
add two more parameters Vtbh"^ and Qmh'^ when we use the BAO data. In [22], they 
used the priors of Qbh^ = 0.02273 ± 0.00061 and Qch"^ = 0.1099 ± 0.0063. 

From the measurement of the radial (line-of-sight) BAO scale in the galaxy power 
spectra, the cosmological parameters were determined from the measured values of 

A.,,o(.) = !^MrA^ (11) 

c 
at two redshifts z = 0.24 and z = 0.43, which are Azbao{z = 0.24) = 0.0407 ± 0.0011 
and AzBAoiz = 0.43) = 0.0442 ± 0.0015, respectively [23]. Therefore, we add x^ with 

,U.(p.f!.^^^) = fA...o(0,24)-0,0407y /A...o(0,43)-0.0442^ 

When we add these BAO data to the fitting, we also need to use the parameters Qbh'^ 
and O^/i^ xhe values fi^/i^ = 0.02273 ± 0.0066 and n„^h'^ = 0.1329 ± 0.0064 were used 
in [23]. 

In addition to the above two BAO data sets, the BAO A parameter [30] is usually 
used. The BAO A parameter is defined as 



.(12) 
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and it was measured to be A = 0.493 ± 0.017 [21j, so we add x with 

2 , , /A-0.439\2 
XBAOaKV) = y QQ^^ j • (14) 

Note that the BAO A parameter depends on the model parameters p only; it does not 
depend on the baryon density Qbh"^ smd the Hubble constant h. Although the radiation 
density depends on h, the contribution to the Hubble parameter E{z) is negligible at 
the redshift z = 0.35, so we can neglect the radiation component when we use the BAO 
A data. 
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Both the SN la and the BAO data measure the distance up to redshit 2; < 2; we 
need to consider the distance at hight redshift in order to determine the property of 
dark energy. Therefore, we implement the WMAP5 data. To use the full WMAP5 
data, we need to add some more parameters which depend on inflationary models, and 
this will limit our ability to constrain dark energy models. So we only use the WMAP5 
measurements of the derived quantities, such as the shift parameter R{z*), the acoustic 
scale Ia{z*) and the decoupling redshift z* , to obtain 

XCMB = AXiCoV2"^(x,, Xj)AXj, (15) 

where the three parameters Xj = {R{z*), Ia{z*), z*), Axj = Xi — xf"^ and Cov2(a;j,a;j) is 
the covariance matrix for the three parameters [21] . The shift parameter R is expressed 
as 

\/rr~ I I — /■■2* dz \ 

i?(z*) = ^^^sinn K/|fifc| / -—-= 1.710 ±0.019. (16) 

The acoustic scale I a is 

Uz'^) = l^^^^jf] = 302.1 ± 0.86, (17) 

(1 + z*)rs[z*) 

and the decoupling redshift z* is fitted by [31] 

z* = 1048[1 + 0.00124:{nbh^y°-^^^][l + gi{n^hy^] = 1090.04 ± 0.93, (18) 

_ 0.0783{Qbh^y°-^^^ 0.560 

^' ~ 1 + 39.5(^/^2)0.763' 92 - l + 21.1(f]b/l2)l-81- ^ ^ 

In [21], it was found that f^fe/i^ = 0.02273 ± 0.00062 and h = 0.719i[J:l]i?. 

The SN la data, the BAO data and the WMAP5 data use the distance measurement 
to determine the cosmological parameters. To get the distance scale, we need to integrate 
the equation of state parameter w{z) twice, so the process of double integration smoothes 
out the variation of equation of state parameter w{z) of dark energy. To alleviate the 
problem, we add the Hubble parameter H{z) data. The Hubble parameter H{z) at 
nine different redshifts was obtained from the differential ages of passively evolving 
galaxies in [25], and three more Hubble parameter data H{z = 0.24) = 76.69 ± 2.32, 
H{z = 0.34) = 83.8 ± 2.96 and H{z = 0.43) = 86.45 ± 3.27 were determined recently in 
[26] . Therefore, we add these H{z) data to x^'- 

x?,(p,/.)=i:i^M^#^^^^, (20) 

where ahi is the la uncertainty in the H{z) data. The model parameters p are 
determined by applying the maximum likelihood method of x^ fit- We use the publicly 
available MINUIT code for minimization and contour calculation [32]. Basically, The 
model parameters are determined by minimizing 

2 2 I 2 I 2 I 2 I 2 I 2 /"OT \ 

X — Xsn ' XSaoa ' XBao2 ' XSaoz ' XcMB ' Xh- \^^) 

For the convenience of numerical fitting, we take Qhh"^ = 0.02273 determined from the 
WMAP5 data [21] • For the Hubble constant h, two different values were observed. 
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The Hubble key project found that h = 0.72 ± 0.08, and recently Riess et al obtained 
h = 0.742 ± 0.036 by using a differential distance ladder method [3lj. To account for 
the uncertainty of the Hubble constant, we treat it as a free parameter and then fix it 
at its best fit value. 



3. ACDM model with curvature 

For the cosmological constant, the equation of state parameter w = p/p = —1, and the 
energy density px is a constant. In a curved ACDM model, the curvature term k ^ 0, 
ordinary pressureless dust matter, radiation and the cosmological constant contribute 
to the total energy. The Friedmann equation is 

E{z) = ^ = [n,{i + zf + nui + z)^ + a(i + z)* + n^f', (22) 

-no 
where the Hubble constant Hq = 100hkms~^Mpc^^, Qm = {^t^G pm) I {^Hq) is the 
current matter component, the current radiation component VL,^ = {^ttG pr) / {^Hq) = 
4.1736 X 10~^h~^ [24], the current curvature component Qk = ~^/(c^o-^o) ^^^ ^a = 
1 — Qm — ^k — ^r- In this model, we have two parameters p = {Qm, ^k) and one 
nuisance parameter h. For the fitting to the SN la data, the contribution to the Hubble 
expansion from the radiation is negligible and we usually neglect the radiation term. 





Figure 1. The la, 2a and 3a joint contour plots of i7,„ and ft^ (ilk) for tlie curved 
ACDM model. The straight line in the left panel denotes the flat ACDM model. 



By fitting the ACDM model to the above observational data, we get x^ = 482.13, 
Qm = 0.280to'oi3 and Qk = 0.001 ±0.005. By fixing the nuisance parameter h at its best 
fit value h = 0.702, we obtain the contours of Q^n and Qk- The joint contour plots of 
Qm and Qk or ^a are shown in figure [H Compared with WMAP5 fitting results [24J , we 
find that the current data make a little improvement on the constraints of Qm and Q^- 
The improvement is due to more SN la and BAO data in addition to the H{z) data. 
The result tells us that the flat ACDM model is consistent with current observational 
data at the la level. 
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Table 1. The marginalized estimates of the model parameters in CPL and JBP 

models. 



Model 


x' 


Qrn 


Qk 


Wo 


Wa 


CPL 


481.27 


U-^'y-0.008 


0.005t°:°?? 


J-'UJ-o.oe 


0.47t?:l 


JBP 


481.46 


n O0 1+0.015 

U-281_o.oi 


o.ooot'oZ 


-0.96t°:l 


-0.6111 



4. CPL parametrization with curvature 

In order to investigate the equation of state of dark energy for a curved universe by 
observational data, in this section we study the popular CPL parametrization [13] 

WaZ 



w{z) = Wo + 



l + z' 



(23) 



The dimensionless Hubble parameter including the contributions from dark energy. 
ordinary pressureless dust matter and radiation is 

H{z) 
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where the dimensionless dark energy density is 

nr)E{z) = {l-n^-nk- fi,) x (l + ^)3(i+"'o+^a) exp[-3waz/{l + z)\25) 

In this model, we have four model parameters p = {Qm, ^k, uio, Wa)- By applying the 
observational data discussed in the previous section to the CPL model, we are able to 
get the observational constraint on the model parameters p = {flm, ^k, wq, Wa). The 
best fit is x^ = 481.64, fi„ = 0.278, Q^ = 0.006, Wq = -1.04, Wa = 0.42 and h = 0.70. 
By fixing the parameters Qm, ^k and h at their best fit values, we obtain the contours of 



Wo and Wa and they are shown in figure 2(a) From figure 2(a) , we see that the ACDM 



model is excluded by the observational data at more than 3a level. As we discussed in 
the previous section, we see that the ACDM model is consistent with the observational 
data. The totally different conclusions suggest that the simple x^ error estimation by 
fixing other parameters at their best fit values has some drawbacks because we neglect 
the correlation effects of the other parameters. The degeneracy between parameters 
was not considered in the above method. When the parameters are strongly correlated, 
the error of some parameters will be under-estimated if we fix the other parameters at 
their best fit values. To verify this point, we apply the MCMC method to constrain the 
parameter space p and the nuisance parameters h and Qbh"^. Our MCMC code [28] is 
based on the publicly available package COSMOMC [35]. By using the MCMC method, 
we get x^ = 481.27; the marginalized la errors are Qm = 0.279tn-nn«, ^k = 0.005+°°°^ 



Wo 



-imt'^om and Wa 



0.5: 



-0.3 
-1.5- 



0.008) ^'A: — >J->J>J^-0.011) 

These results are summarized in table [H The 



marginalized la, 2a and 3a contour plots of Wq and Wa are shown in figure 2(b) From 



figure 2(a) and 2(b) , we find that the marginalized la contour obtained by using the 



MCMC method includes the 3a contour in figure 2(a) From figure 2(b), we see that 



the ACDM model is consistent with the CPL model at la level. 
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Figure 2. The la, 2cr and 3t7 contour plots ol wq and Wa for the eurved CPL 
parametrizatfon. '+' denotes the point corresponding to the ACDM model, (a) Joint 
contours of Wq and Wa by fixing the other parameters at their best fit values, (b) 
Marginalized contours of wq and Wa obtained from the MCMC method. 





0.26 0.28 0.3 0.32 



-0.01 0.01 0.02 0.03 





(b) 



Figure 3. The marginalized distributions of the CPL model parameters p are shown 
in (a). The solid lines are marginalized probabilities and the dotted lines are mean 
likelihoods, (b) Reconstructed evolution of w(z) and the shaded areas are la, 2a and 
3cr errors. 



The marginalized distributions of the model parameters are shown in figure 3 (a 



The solid lines are marginalized probabilities and the dotted lines represent mean 



likelihoods. From figure 3(a), we find that the likelihood of fi^ has a local maximum 
around fi^ ~ 0.02. Even we take fi^ = 0.02, the value of x^ is not far from the minimum 
value of x^- Due to the degeneracy between model parameters, if we fix the other model 
parameters at their best fit values, then the joined contours of wq and Wa are under- 
estimated, and the conclusion drawn from the under-estimated contours is not reliable. 



The results in figure 2(b) and 3(a) verify this point. By using the marginalized contours 



of Wq and Wa-, we reconstruct the evolution of wi^z) in figure 3(b), From figure 3(b), we 
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find that w{z) < at more than 3a confidence level up to redshift z = 2, and the ACDM 
model is consistent with the CPL model at the la level. To account for the correlations 
between model parameters, we need to use the marginalized probability. To see whether 
this happens only for the CPL model, we analyze the JBP model in the next section. 



5. JBP parametrization with curvature 

In this section, we consider the JBP parametrization [27] for dark energy with the 
equation of state in the form below 

w{z) = Wo+ 2 ^2 ■ (26) 

(1 + zy 

The corresponding dimensionless dark energy density is then 



^de{z) 



[l-Vtm-^k- fir) X (1 + Z 



,3(l+«>o) 



exp 



^Waz" /2{l + zf (27) 



Wn 



We first use the 



In this model, we also have four parameters p = (fim, fifc? Wq, 
simple x^ method to fit the model. The best fit is x^ = 481.84, fi„ = 0.281, Vtk = 0.0015, 
Wq = —0.97, Wa = —0.03 and h = 0.70. By fixing the parameters Qm, ^k and h at their 
best fit values, we obtain the contours of Wq and Wa and they are shown by the solid 



lines in figure 4(a) Unlike the CPL model, the ACDM model is consistent with the 
JBP model at the la level. To verify this conclusion, we also apply the MCMC method 
to the JBP model. By using the MCMC method, we get x^ = 481.46; the marginalized 



la errors are Qr 



0.281 



+0.015 
-0.01 ) 



n, 



0.0001°:°°^, ^0 



OMZlil and Wa = -0.6; 



-1.9 
-1.6- 



These results are summarized in tabled! The marginalized la, 2a and 3a contour plots 



of Wo and Wa are shown in figure 4(b) From figure 4(a) and 4(b), we see that the 



contours of Wq and Wa are consistent although the constraints from the MCMC method 
are a little larger because we consider the correlations among all the parameters in the 
MCMC method. The ACDM model is also consistent with the JBP model at the la 
level. 

The marginalized distributions of the model parameters p are shown in figure 
The solid lines are marginalized probabilities and the dotted lines represent 



5(a; 



mean likelihoods. From figure 5(a) we find that the probability distributions of the 
parameters are more or less Gaussian. By using the marginalized contours of wq and 



Wa, we reconstruct the evolution of w{z) in figure 5(b), From figure 5(b), we find that 
w{z) < —0.2 at more than 3a confidence level up to redshift z = 2, and the ACDM 
model is consistent with the JBP model at the la level. 



6. Conclusions 

Applying the simple x^ method, we fitted the CPL and JBP models to the combined SN 
la, BAO, WMAP5 and H[z) data, and obtained the constraint on the property of dark 
energy. In both CPL and JBP models, there are four parameters p = (fim, fifc, wq, Wa). 
When we apply the BAO and WMAP5 data, we need to add two more parameters 
Qfyh"^ and h. If we make the joint error analysis, we have six parameters and it will 
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Figure 4. The la, 2a and 3a contour plots of wq and Wa for the curved JBP 
parametrizatfon. '+' denotes the pomt corresponding to the ACDM model, (a) Joint 
contours of wq and Wa by fixing the other parameters at their best fit values, (b) 
Marginalized contours of wq and Wa obtained from the MCMC method. 





0.26 0.28 0.3 0.32 



-0.01 0.01 0.02 





(b) 



Figure 5. The marginalized distributions of the JBP model parameters p are shown 
in (a). The solid lines are marginalized probabilities and the dotted lines are mean 
likelihoods, (b) Reconstructed evolution of w(z) and the shaded areas are la, 2a and 
3(T errors. 



be hard to get a good joint constraint on all these parameters. Therefore, we take 
flbh"^ = 0.02273, and find out the best fit values of the parameters p and h which 
minimize x^; then we fix the parameters flm, ^k and h at their best fit values to obtain 
the joint constraints on wq and Wa- For the CPL model, the contours of Wq and Wa (see 



figure 2(a)) show that the ACDM model is excluded at more than 3cr level. The JBP 
model is consistent with the ACDM model at the la level. Since we get the contours 
of Wq and Wa by fixing the other parameters at their best fit values, we neglect the 
correlation effects of the parameters and the conclusion based on this method may not 
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be reliable. To confirm this, we use the MCMC method to analyze the CPL and JBP 
models and obtain the marginalized probabilities of the parameters. For the CPL model, 
the probability distributions of Qk and Wa are skew distributions, and the marginalized 
la errors are n^ = 0.279to:ooi, ^k = 0.005+^??, Wq = -1.05t[]:g^ and Wa = 0.5i?i In 
the CPL model, the probability distributions of ^2^ has a local maximum in addition 
to a global maximum. The uncertainties in Qk and the degeneracies between fik, wq 
and Wa lead to under-estimation of the error contours of Wq and Wa if we fix Q^ at its 
global best fit value, and the wrong conclusion that the ACDM model is excluded at 
more than 3a level. However, this does not happen for the JBP model. For the JBP 
model, the parameters have Gaussian distributions, and the marginalized la errors are 
n^ = 0.281+°:°^, fifc = 0.000t[J:°|]^, Wo = -OMttil and Wa = -o.Qtll 

In summary, in addition to use the usual SN la, BAO A or BAO distance ratio, 
and WMAP data, we also use the radial BAO measurements and the H{z) data to fit 
the CPL and JBP models. We find that the equation of state parameter of dark energy 
w{z) < at more than 3a level in the redshift range < z < 2, and the fiat ACDM 
model is consistent with the current observational data at the la level. Furthermore, 
we find that we need to do the marginalized analysis to estimate the errors of the model 
parameters. 
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